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Abstract — An advanced transitional turbulence model 

( ^ ) has been applied to the simulation of the typical 

hypersonic flows in the present investigation, which is corrected 
with empirical correlation for the transition consideration. 
Combining PSE (stability analysis equation) with eN prediction 
method, it has been studied how to find out the position of the 
most unstable flow within the boundary layer, so that the 
transition empirical correlation parameters would be more 
appropriate for the specified flow while being independent of 
designer’s experiences. The typical hypersonic flow field over a 
blunt cone has been simulated, through which it could be found 
that both the computational results of the wall heat flux and 
boundary layer transition region agree well with the 
experimental results. Besides t, both the interaction between 
shock wave with the boundary layer on a flat plat, and the flow 
across a compression corner have been simulated by the same 
model as well. It was found that the present methods are also 
applicable for the simulation of complicated flow with shock 
wave/boundary layer interaction, and the results turned out to 
be superior to those employing traditional turbulence models. 

Index Terms — Hypersonic flow, PSE Stability Equation, 
Transition, Shock wave, Boundary layer 


i. Introduction 

Turbulence flow is an important phenomenon in fluid 
dynamics. In the aerodynamic design one of the most 
difficulties is to simulate the complex turbulence flow. Direct 
numerical simulation (DNS) of turbulence and Large Eddy 
Simulation (LES) have good theoretical basis and broad 
application values. However, limited by the current 
state-of-the-art calculation conditions, they are not applicable 
to engineering modeling yet. Though DNS is the most 
accurate and appropriate method for transition and turbulence 
study, it still strongly constrained by computer resources and 
algorithmic limitations. Most current DNS studies are only 
conducted for very simple geometries. Therefore, using the 
current turbulence models to solve engineering turbulence 
problems is still the most effective and most practical choice. 
Hypersonic flows involve complex flow phenomena, such as 
large adverse pressure gradient, flow separation and 
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re-attached as well as the strong effect of compressibility. In 
the flow transition region, aircraft's surface heat flow rate and 
the surface friction coefficient will increase dramatically. 
Therefore, to accurately predict the location of the transition 
starting point and the length of transition zone is engineering 
is the key point of engineering design, especially in the 
thermal protection system design. When calculated by 
turbulence model, it is obviously unreasonable if the 
transition location is directly given by empirical 
formula[l,2]. It is generally believed that the transition and 
instability mechanism of three-dimensional hypersonic flow 
have three ways: second-mode, cross-flow and bypass[3]. In 
order to study these mechanisms, researchers have been 

carried out a great deal of research work. The ^ model 
proposed by Langtry and Menter[4,5] overcomes this 
difficulty by predicting transition only by the use of local 
quantities. Martine had attempted to investigate a well-suited 

hypersonic intake using innovative ^ model in 

combination with the in-house correlations for transition 
onset and length. The new model has a great improvement 
than those used before but it still needs to determine the 
Reynolds Number of transition starting point by empirical 
correlation and increase the uncertainty of the numerical 
calculation [6] . 

Stability analysis is one means of studying the physics of 
the transition problem, such as the LST and PSE methods [6] . 
The most popular theoretical prediction method in practical 
use is to calculate the development of disturbances, and give 
transition position to reveal the mechanism of it, which is 
based on linear stability theory of the e N prediction method. 
For example, a Hypersonic Boundary Layer package has been 
developed for the stability and transition analysis [710] . It 
generates a prediction for the evolution of an initial 
disturbance when marching downstream from starting point 
to predict the onset of transition, an experimental correlation 
is required through the semi-empirical e N model However, its 
limitations is so large that it could not be completely used in 
the engineering designs [1U2] . It is not so clear about these 
mechanisms because of the complexity of the hypersonic 
flow. The numerical results and experimental ones show a 
large deviation. 

The paper presents an approach to joint these two methods. 
The specific way is to use e N method to estimate boundary 
layer instability location, to calculate the momentum 
thickness Reynolds number of the transition starting point, to 
amend transition model empirical relations according to the 
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situation, to provide the experience factors that could be 
directly used to flow field calculation and to reduce reliance 
on the designer experience. The results show that this method 
can more accurately simulate the transition process and the 
shock wave/boundary layer interaction caused by separated 
flow. The wall pressure and heat flux of the transition zone 
and separation zone are in good agreement with the 
experiment, which gives a more accurate simulation of 
supersonic or hypersonic flow of the complex multi-wave 
system, transitional flows, and aerodynamic heating and so 
on. 


H. r Re « TRANSITION MODEL 

To accurately solve hypersonic flows, the time-averaged 
Navier Stokes equations are solved by spatial discretization. 
Fluent software is adopted as the numerical solver. MUSCL 
scheme with third-order accuracy and implicit time-advance 
scheme are also used. 

The flow field changes severely in the shock 
wave/boundary layer interaction region, where the turbulence 
is unbalanced, especially under the interaction of strong 
shock wave. However, turbulence model treats the turbulent 
fluctuation statistics and the averaged quantities separately, 
and uses the averaged quantities to model the pulse 
statistics [13]. The general Turbulence models are not enough 
to accurately reflect the characteristics of the 

non-equilibrium, that is said, have no adequate abilities to 
describe the dependence of the pulse statistics quantities on 
the averaged in the process from the non-equilibrium to the 
balance. Thus the physical characteristics of turbulence 
pressure enhanced after the shock wave as well as the heat 
flow are all overestimated[14]. In this paper, the improved 

^ 6 two-equation turbulence model as well as the UDF 

function of Fluent are both adopted to control the 
experimental correlation formula in the transition model and 
better solutions of this problem are achieved. 

Based on Menter's SST model, Langtry and Menter[4-6] 

y — 

put forward the 7 6 model which could associate the 

transition starting point and the length of transition zone. Two 
transport equations about transition intermittent factor and 
transition initiation thickness Reynolds number were 
introduced. 

(1) Intermittent factor transport equation: 


^ + £(P^y) = _a_ [(|i+ 1 L) _a L] 

dt 8xj 7 7 8xj a f 8xj 


( 1 ) 


(2) Transition initiation thickness Reynolds number 
transport equation: 


d(pRe e ,) g(pt/ y Re e ,) n 

dt + 8xj 9t 


d dRe fi , 

+ — K,(p + p,)— -] 

OXj OXj 


( 2 ) 


Among the source terms of the transition, ^ e ° c 
is used to associate with the starting point and 
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^length /(^ e er) - s usec j tQ assoc j a t e w ith the transition 
Re 

length. 01 is the transition initial Reynolds number and is 

F 

given by experiment relation formula in Fluent solver. 0/ is 
blending function to turn off the boundary layer source term 

Re 

and the control of transport options of e U 

The solution is firstly performed with conventional 

^ - Re 6 turbulence model and the initial data of flow field are 
obtained for the following PSE solution. Then the gained 
outer boundary layer parameters are passed to the PSE solver 
for boundary layer flow stability analysis. In this step, the 
most critical issue is to use the location of the starting 
instability point and the transition point gained from 
boundary layer analysis to modify the critical transition 

Reynolds number 0C and the initial transition Reynolds 

number Qt based on the formulate(3)-(6) in reference [13]. 
This step is used to replace the experience coefficients of 
Fluent transition model to close the transition model. Finally, 
using UDF function, substituting the revised transition 
formula into the Fluent transition factor and momentum 
thickness Reynolds number of the transport equation, solving 
the flow field until it convergences as is shown in figure 1. 


Re flc = 


^ 

(-0.0427m 3 + 0.42337w 2 + 0.01 1 8Tu x + 1 .0744) ' (3) 


Re 0 , = [1 173.51 - 589.4287w + 


0.2196 


Tic 


]■ F(k e ) ,Tu<l.3. 


Re fl , =331.5-[7’m-0.5658] -0 ' 671 ■ F(T b ), 7m >1.3. 

(4) 

Tu 

F(\) = l - [— 12.986X 0 - 123.667./ - 405.689X e 3 ]-e 05 , \ > 0. (5) 

t inn >/2*73 pQ 2 dU 

U ju ds 


( 6 ) 


Which are illustrated in detail in referenced 3]. While 
coupled with Fluent, a code for UDF is compiled in the 
following form: 

#include "udf.h" 

DEFINE_TRANS_FLENGTH(user_Flength, c, t) 

{ 

...// 

} 


DEFINE_TRANS_RETHETA_C(user_Re_thetac, c, t) 

{ 

....// 

} 


DEFINE_TRANS_RETHETA_T(user_Re_thetat, c, t) 

{ 

...// 


Reference [13] used these relation equations to solve 


Re e c 


and F ‘ ength 


in order to close transition model equations. 
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However in this paper, PSE boundary layer stability analysis 
method is applied to deal with them. 

m. PSE STABILITY EQUATIONS 1151 

Before applying the nonlinear boundary layer stability 
analysis to the methodology developed in this paper, the open 
forms of Linear Stability Theory (LST) and Parabolized 
Stability Equation (PSE) are introduced firstly. To derive the 
LST, the flow is assumed to be locally parallel. Substituting 
the disturbance quantities into the linearized Navier Stokes 
equations, it can be obtained the first order differential 
equations of LST in matrix form: 


A +v ^ + v„A +v A t v,' = 

dt as an ar 

I *-\2 I *-\2 I o 2 I *-\2 I o 2 I 

v 8 -± + v d -± + y d -± + y e -± + y 8 -± + y 8 * 

J cc /% 1 nn /% ’ rr /% ’ T cr T nt 


(14) 


a* 


dn 


dr 


ds dn 


dsdr 


dndr 


In which ^ ^ and i} are all 5x5 operator matrixes. 

Based on the assumptions that the disturbances are 
composed of a fast-oscillatory wave part and a 
slowly-varying shape function, the disturbance equation 
above is simplified into the following form: 


q (Xp^ ,Xyt) = qe 


i(ax j + (tit ) 


dq , 


</(/| du^ 


dT 


q = [u 1 ,- — ,m 9 ,- — ,T,~ — ,u?,p] . 
1 dxn, dx~ dx~ D 


(7) 

( 8 ) 

(9) 


d d 

+> 77 + A ~ + A) 

dn dn 


. B di|> v ss ^ v sn 

t), = 21 — (15) 

Re dx Q Re dx c Re dxdn 


This is an eigenvalue problem and can be solved through 
numerical discretization method. The boundary conditions 
are as follows: 


In general the three dimensional boundary layer are not 
parallel, substituting the disturbance forms into the linearized 
Navier-Stokes equations, it can be obtained the second order 
differential equations of PSE written in matrix form: 


q (x l ,x 2 ,x 3 ,t) = q{x l ,x 2 )e 


i{a(xi +(3x3 +(D 0 


A^- + A ^ 

2 


d 2 q 


dx 2 2 


dx x x 2 dxy 


C l ^- + C 2 ^- + I\q = 0 


' dx 0 


q - [p,w 1 ,w 2 ^3 5 T ] 7 


( 10 ) 

(ID 

( 12 ) 


Then the nonlinear PSE equations are developed based on 
the assumption of small perturbation as illustrated in this 
paper. Nonlinear PSE equations are based on the assumption 
of small perturbation and can be directly used to solve 
three-dimensional problems. The fluctuations of pressure, 
density, temperature, viscosity, thermal conductivity and 
other physical quantities are taken into account in 
calculation [7]. 

In this paper, the disturbance variation in curved 
coordinate is adopted as follows: 


q =q(s,n)ex p i 


J ^a(s)ds + fi' 


z + (ti-t 


(13) 


u =v = w = f = 0 - y = 0 ; ( 16 > 

u = v = w = f = 0 ’ y— >oo. (17) 

Since the disturbance evolutions are mainly focused on the 
stream- wise and normal direction for specific problems, the 
PSE analysis system above is also suitable for prediction of 
two-dimensional or axi-symmetric flows. For the 2D and 
axi-symmetric flow models discussed in the following, a 
finite thickness is assumed in the cross (or the span-wise) 
direction, in which the disturbance quantity is set to be 0 
during the 3D solver marching process. 

In order to accurately solve the boundary layer equations, 
sufficient number of grid points is needed within the 
boundary layer calculation region. In this paper, 200-300 
points can be arranged in the normal direction of the wall and 
high-precision spectroscopy method is used to solve the 
discrete boundary layer disturbance equations. During the 
solving process of the nonlinear PSE equations, the solution 

variables and ^ needs to be updated. 



In which, ^ [p > u > v > w » ^ ] [ s disturbance 

quantity. & and @ are respectively the wave number of flow 
direction and radial direction. &) i s the disturbance frequency 

and s is the curve coordinate along flow direction. ^ S,n ^ is 
the disturbance amplitude. This perturbation is substituted 
into the Navier Stokes equations and the mean flow is 
subtracted, then the nonlinear disturbance equation in 
orthogonal curvilinear coordinates is as follows: 


The introduction of energy conservation condition on the 
above equation can eliminate the correlation of n: 


E = J^p(|w'| 2 +|v'| 2 +|w'| 2 )d?z 


(19) 


The superscript "+" in the equation indicates complex 
conjugate. So the transition prediction can be done by 
accumulating disturbance growth rate along flow direction. 
Its integral form is as follows: 


N( co, s) = f ' [-Im(a) + — —tf s 
Js o 2 E ds 


( 20 ) 
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In which, is the first position of a neutral point within a 
given frequency. In this way, the occurrence of transition can 
be predicted by the size of N values. 

A detailed derivation of the above equations and the 
validation of PSE calculation process can be seen in 
Reference [15]. 

IV. RESULTS AND DISCUSSION 
A. SE Application 


Outer Boundary Layer Parameters 



Fig.l. Schematic diagram of the solution process 


In this paper, the solution schedule mainly includes the 
following three steps .First, the model is preliminary 
calculated using Fluent software and the transition models 

(for example the r 6 model) using conventional common 
empirical coefficient. Then, the gained outer boundary layer 
parameters are passed to the PSE solver for boundary layer 
flow stability analysis. In this step, the most critical issue is to 
use the location of the starting instability point and the 
transition point gained from boundary layer analysis to 

Re 

modify the critical transition Reynolds number 0C and 

Re 

transition initial Reynolds number 0r based on the 
formulate(3)-(5) in reference [13]. This step is used to replace 
the experience coefficients of Fluent transition model to close 
the transition model. Finally, using UDF function, 
substituting the the revised transition formula into the Fluent 
transition factor and momentum thickness Reynolds number 
of the transport equation, solving the flow field until it 
convergences. Fig.l illustrate the solving process. 

The grid sensitivity for the numerical simulation of 
transitional boundary layer has been studied in detail[14,15], 
which point out that the Menter & Langtry(SST Gama) 
turbulence model has better ability to the simulation of 
boundary layer transition. Through the specification of 
transitional onset momentum Reynolds number and inlet 
viscosity ratio from test data or designer experience, it is able 
to provide better simulation of heat transfer within 
transitional boundary layer. Specially the SST model has 
higher requirement for the numerical grid density with the 
first layer grid nodes near wall surface satisfying y+<l . In this 
paper all cases are computed on the grid with 0.05<y+<0.5 for 
the first grid layer near the wall surface. While for the grid 
sensitivity for the PSE solver, it is performed during the TS 
wave stability analysis for the specified simulation, where 


only the absolutely similar wave numbers obtained on two 
different grids are considered to be confidential. 

The current studies [2] have shown that when the incoming 
flow Mach number is less than 4, the flow instability of the 
focus mode is a first-order mode (TS mode). When the Mach 
number is greater than 4, the master transition process is the 
second-order mode disturbance. For the cone hypersonic flow 
model, due to a relatively low velocity after bow shock, it is 
difficult to find the instability second-order mode. Therefore, 
to study the instability of boundary layer flows, this paper 
introduce the PSE equations to calculate the first-order mode 
wave (TS wave) while the initial disturbance frequency is 
200khz. For the application of TS wave to the transitional 
hypersonic flow in this paper, it is attempted to predict the 
most unstable location of the boundary layer, the precise 
disturbance amplification rate maybe not necessary. Only the 
approximate unstable locations should be specified to 
provide the key transitional parameters such as the critical 

transition Reynolds number 0c and the initial Reynolds 
Re 

number Qt , which are needed to complete the coupling 
solution with mean flow code. So only TS wave( the 1st 
mode) is considered for simplification and the conclusion is 
validated. 

B. Case 1[16] 



Fig. 2. Geometry Model of a Cone 


Initial mean flow solution 

Fig. 2 is the geometry model of a cone, The cone angel of 
the model is 7 degree, cone length is 2.35m and the head 

radius is 5mm. The inflow Mach number ^°°is 7.15, inlet 

p T 

pressure 00 is 7722Pa and inlet temperature 00 is 214K. 

T 

The temperature at cone wall ^is 300K. Two-dimensional 
axisymmetric calculation method is applied. y+ is less than 

0.5 and ^ Uqo is about 0.01%. 

The solution is firstly performed with conventional 

y ~ Re 6 turbulence model and the initial data of flow field are 
obtained for the following PSE solution. 

Eigen-value solution on two sets of grid configurations 
Fig3 displays the real and imaginary part of stream-wise 
wave number a at different axial positions of 0.2 and 0.85. 
As is shown in the figure that the eigenvalue solution is 
performed on two sets of grid configurations, with the grid 
nodes 300 x 151 and 300 x 251 respectively. Let nx represents 
the grid node number in streamwise direction(parallel to the 
wall), and ny in normal direction(perpendicularly to the wall). 
It is considered in this paper that the same eigenvlaue solved 
under normal grid nodes of ny=151 and ny=251 is reasonable 
for PSE solution, and the eigenvalue with the minimum 
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imaginary part represents for the one which leads to the most 
instability of the flow. 

The current studies [2] have shown that when the incoming 
flow Mach number is less than 4, the flow instability of the 
focus mode is a first-order mode (TS mode). When the Mach 
number is greater than 4, the master transition process is the 
second-order mode disturbance. For the cone hypersonic flow 
model, due to a relatively low velocity after bow shock, it is 
difficult to find the instability second-order mode. Therefore, 
to study the instability of boundary layer flows, this paper 
introduce the PSE equations to calculate the first-order mode 
wave (TS wave) while the initial disturbance frequency is 
200khz. For the application of TS wave to the transitional 
hypersonic flow in this paper, it is attempted to predict the 
most unstable location of the boundary layer, the precise 
disturbance amplification rate maybe not necessary. Only the 
approximate unstable locations should be specified to 
provide the key transitional parameters such as the critical 

Re 

transition Reynolds number 0C and initial Reynolds 
Re 

number 0/ , which is needed to complete the coupling 
solution with mean flow code. So only TS wave( the 1st 
mode) is considered for simplification and the conclusion is 
validated. 


Mean flow solution coupled with PSE analysis results 
Fig.4(a) shows the wave-number distribution along flow 
direction. The starting location of the instability point is 
located at x=0.35, the location of the most unstable point is at 
x=0.65 where the N value is 7.2 here, as is shown in Fig.4(b). 
According to the formula (3) the calculation the critical 

Re 

transition Reynolds number 0C and initial Reynolds 
Re 

number Qt on these two positions is carried out and finally 

determines length . So that these three parameters can be 
re-substituted into transport equations, through the UDF 
compile processed, continue to iteratively solve the whole 
flow field. The wall heat flux distribution is shown in Fig. 5. 
As we can see, there is a serious deficiency in transition 
simulation process while applying traditional method. If 
transition model is used directly, the length of the transition 
zone, as well as the wall heat flux of the transition process 
have been predicted a larger deviation. The transition zone is 
too small and the heat flow is significantly high. The 
empirical parameters of the transition model need to be 
constantly adjusted. When using this transition amendment 
method of this article, because it can more accurately give the 
starting transition Reynolds number, the occurrence and 
development of the transition process simulated are more 
rationalized. The lengths of the transition zone and heat flux 
after transition are more close to experimental results. This is 
very useful to simulate a complex simulation of hypersonic 
wave structures, separation of boundary layer flows and 
aerodynamic heating phenomena. 

It should be pointed out that the value of N used in this 
article only was an accurately one as long as possible. That is 
said, if the obtained transition momentum thickness Reynolds 
number is more accurate, the calculation of the transition 
process can be improved. This greatly reduces the designers' 
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dependence on empirical parameters. Currently one of the 
problems is the PSE analysis process requires a different 
density level of the grid, through the shooting method to 
derive characteristic matrix eigenvalues, to find the minimum 
and the imaginary part of eigenvalues at different grids 
overlap levels as an expression of flow field in the most 
unstable point. Thus, when the grid point increase n times, the 
matrix size will be n x n times and the computing time will 
also increase n x n times. The PSE actual time-consuming of 
this article roughly is equivalent to a flow field calculation 
time. How to improve the speed of PSE solution will be 
discussed in the next step. 



-3- 


-5' 1 1 1 1 1 ' 1 1 1 ' 

-1 -0.5 0 0.5 1 1.5 2 2.5 3 3.5 4 

Real part 

(a) s=0.2 


t * ■ 

* 
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* 

* * * 

* * * 
* 

* 


* 


x ny=12l 
+ ny=251 


1.5 2 

Real part 


(b)s=0.85 

Fig. 3. The real and imaginary part of stream- wise wave 
number a at different axial positions(a) s=0.2, (b)s=0.85 
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Fig.4 . Distribution of stability analysis factors on Cone Fig. 6 . Numerical Calculation Model 

Surface (a)Wave Number 



Fig.4. Distribution of stability analysis factors on Cone 
Surface, (b) Integrated N Vales 


c Case 2 1171 

Fig . 6 is the numerical calculation model. The supersonic 
flow passes through an inclined board, creates a shock wave 
by compression effect. The shock incidents to the plate 
boundary layer and thus leads to a shock wave/boundary layer 
interference. The ramp length is 8.89xl0-2m and the flat 
length is 30.48x1 0-2m. The plane angel a and leading edge 
coordinate C are shown in Table. 1. The inlet conditions are as 
follows: 

R e L = l-5xl0 5 , M aoo = 2.96 , Tu„ =0.2%, 

T = 247/T anc j P ~ 4-07 x 10 the wa |] s are treated as 
adiabatic. The grid nodes for mean flow calculation are 

200 x 200 , and y+ <l . 



Ma = 2.96 


& 


yc) 


~TT / / A 

y 


//// / ' ///?// A / 7-X 


Tablet. Plane Angel a and Leading Edge Coordinate 


(* c ,y c ) 

a/(°) 

x c /10 2 m 

y c / 10 2 m 

9.39 

-9.144 

4.674 

12.27 

-9.500 

5.436 


As is shown in Fig. 7 the wall pressure numerical results are 
compared with the experimental values. The pressure 
suddenly rises through the shock wave. After the shock wave, 
the subsonic flow passes back through the boundary layer 
against the mainstream so that the pressure upper of the shock 
incident wave rise. So boundary layer separation occurs 

* * * • • » r \f Pp 

before the shock incident point. At this situation, r 0 
transition model of Fluent Software is used to accurately 
obtain the value of pressure increase after shock wave. 
However, the location of separation point and re-attach point 
and the length of separation zone are still having a certain 
degree of deviation. Especially for the situation when a 
equals 12.27degree, the predicted separation zone is a little 
lager and the separation point is obviously moved upwards. 
By using the methods above, the boundary layer stability 
analysis is carried out while N is selected as 8.0. The initial 
Reynolds number at the boundary layer instability transition 

point is gained. Applying the revised r 6 transition 
model, an accurate pressure rise is obtained ( oc = 9.39 , 

pfPo~ 3 35 w hile a = 12.27 , P^Po~ 4 9 y The 

prediction of the location of separation point and re-attach 
point and the length of separation zone has a good fit with the 
experimental values. In particular, the location of the pressure 
sudden increase, the pressure of the pressure platform area 
and the subsequent pressure stable value are well agreed with 
experiments. It indicates that the solution of this article is 
better than those using traditional turbulence model. See from 
Fig . 8 and Fig. 9, the calculation method of this paper 
accurately calculates the shock wave structure that contains 
shock incident to the laminar boundary layer on a plate after 
the shock wave reflection, expansion, and re-compression. 
The calculation is in good agreement with the experimental 
results. 

Figure 10 shows the distribution of the axial component u 
along the normal direction of the wall, in which xl and x7 
denote the section surface outside of the separation region, x 2 
and x 6 represent section surfaces passing through the 
separation and reattachment points, and x3, x4 and x5 stand 
for the sections within separation region. From the figure it 
can be seen the significant influence of the shock wave on the 
velocity distribution of boundary layer. 

D.Case 3" 81 

Hypersonic axial symmetric compression corner flow with 
a skirt section is a typical example to verify the CFD 
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prediction capacity of hypersonic flow surface pressure and 

T =31 IK 

heat flux. The inlet flow conditions are: w , 

=7.05 f Px = 516.0 pa and T x =81.2 K The skew 

angel of the skirt section is 35 degree. The thickness of 
boundary layer is 2cm before the corner separation point. 
Two-dimensional axial symmetric calculation is adopted. 

The grid nodes are 200 x 200 and < ^ . 

Fig. 11 is a typical geometry model of corner compression 
in hypersonic flow. N is selected as 8.0 when calculating 
boundary layer stability by PSE method. Joint with flow field 
solution, wall pressure has been shown in Fig. 12 and wall 
heat flux density shown in Fig. 13. In these figures, the results 
are compared with the experimental values and those do not 
use PSE for transition model updating. By comparison you 
can see, the application of PSE on amendment of transition 
gets a more accurate pressure distribution and heat flow rate 
peak, especially for the peak position forecast. At the same 
time, the heat flow overshoot phenomenon right after the 
peak pressure is described more clearly. It has obvious 
advantages in comparison with no amendments to the 
transition model. This is very meaningful to the heat 
protection design of hypersonic aircrafts. 



Fig. 7. Pressure distribution on Flat Surface (a) a =9.39° 



Fig. 7. Pressure distribution on Flat Surface (b) a =12.27° 




Fig. 9. Counters of Density in Interaction Zone (Numerical) 



Fig. 10. Distribution of axial velocity in the interaction zone 



Fig.l 1. Geometry Model of Corner Compression in 
Hypersonic Flow 
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Fig. 14 shows the density and pressure isolines within the 
interaction region, in which a represents a weaker shock wave 
caused by the concave of the stream lines near the separation 
point of the boundary layer, b is a strong shock wave from the 
reattachment point in separation region, which interacts with 
shock a, forming a stronger shock wave c. The flow behind 
shock wave b comes from the free incoming flow only 
compressed by shock wave c, so the flow energy loss behind 
shock wave b is not large, with higher local pressure there. 
Therefore a set of expanding waves are created, leading the 
pressures behind shock wave b and c to be equilibrium. The 
above illustration of wave structure can be verified through 
the comparison with reference [19]. In the fig. 15 of the 
pressure isolines, it is shown that a higher pressure region 
exists behind the reattachment point, which is corresponded 
to the wall pressure overshoot. The reattachment point is at 
x=0.025m and the pressure overshoot is at x=0.035, the two 
positions are very close. 



Fig. 14. Density distribution in interaction zone 



Fig. 15. Pressure distribution in interaction zone 
V. CONCLUSIONS 

Following the hypersonic flow transition phenomenon, 
the PSE disturbance equations and N-S equations have been 
solved jointly and semi-empirical e N method has been 
employed in the present numerical investigation. The 
objective is to identify the unstable point within the boundary 
layer so that the y - Re 0 turbulence model could be amended 
in Fluent software, correlating the transition relation formulas 
with Reynolds number based on initial momentum thickness, 
the critical Reynolds number as well as the transition zone 
length. It successfully reduced the reliance upon the 
experiments and experiences during the transition region in 
the hypersonic flow process. For the complex transition flow 
around cone surface within the hypersonic flow regime, and 
the interaction between shock wave and boundary layer, the 
capability of y-Re e turbulence model and associated 
guidelines have been proven yielding to more accurate 
supersonic and hypersonic numerical results. 

NOMENCLATURE 

T : Temperature( K) 

/l : Thermal conductivity(W/mK) 
p : Density(Kg/m 3 ) 
p : Pressure(Pa) 

h : Convective heat transfer coefficient(W/m2K) 

U : Flow quantities 
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q,(Q \ Turbulence quantities 
U : Disturbance Amplitude 

U : Mean flow quantities 
S : Stream- wise direction 
U ' : Disturbance quantities 
fl : Normal direction 
CC : Stream- wise wave number 
j3 : Span- wise wave number 
/ : Disturbance frequency 
t : time(s) 

Aq,A l ,A 2 : Matrix 
i : Imaginary part 

Tu ^ : incoming turbulence intensity 

Tu : local turbulence intensity 
y : Intermittent factor 

Re^ : transported momentum thickness Reynolds 
number 

Re^ : initial transition Reynolds number 
Re 0c : critical transition Reynolds number 
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